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Using resonance ultrasound spectroscopy, 
we measured the monocrystal elastic- 
stiffness coefficients, the Voigt C { : , of 
TiB 2 . With hexagonal symmetry, TiB 2 
exhibits five independent Cg : C u , C 33 , 
C44 , C 12 , C 13 . Using Voigt-Reuss-Hill 
averaging, we converted these monocrystal 
values to quasiisotropic (polycrystal) 
elastic stiffnesses. Briefly, we comment 
on effects of voids. From the Cy , we 
calculated the Debye characteristic 
temperature, the Griineisen parameter, 
and various sound velocities. Our study 
resolves the enormous differences between 
two previous reports of TiB 2 's Cy . 
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1. Introduction 

Currently, we see a renascence of research on transi- 
tion-metal diborides. For example, well-known are 
many current studies on ReB 2 and OsB 2 as superhard 
materials, rivaling even diamond. 

Titanium diboride's physical and mechanical proper- 
ties received review elsewhere [1, 2]. This compound is 
known well for low mass density, high hardness, high 
melting point, low electrical resistivity, good thermal 
conductivity, and good chemical inertness. However, 
the problems in preparing full-dense, high-quality 
monocrystals preclude their extensive study. Akimitsu 
and colleagues [3] reported superconductivity at 40 K 
in a same-crystal-structure companion compound: 
MgB 2 . 

The importance of elastic-stiffness coefficients for 
both science and technology also received review [4]. 



Titanium diboride's monocrystal elastic stiffnesses 
appeared in two reports. The values of Gilman and 
Roberts [5] depart strongly from the more recent report 
of Spoor and colleagues [6], mainly in the off-diagonal 
C 12 and C 13 . The first measurements were made by a 
pulse-echo method, the second by resonance ultrasonic 
spectroscopy, but by a very different mechanical setup 
than used in the present study. 

Here, we report a third set of measurements, 
which agree closely with the Spoor et al. results. From 
our Qj measurements, we estimate the Debye charac- 
teristic temperature, the quintessential harmonic-lattice 
property. Also, we estimate the Griineisen parameter, 
the quintessential anharmonic property. 

Figure 1 shows the TiB 2 crystal structure, obviously 
a layered structure suggesting strong elastic anisotropy, 
not observed as shown below, creating a conundrum. 
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Fig. 1. Schematic crystal structure of titanium diboride (A1B 2 type, 
C32, hexagonal, P6/mmm, hP3, M=l, D 6h \ No. 191). Small 
spheres represent boron atoms, large spheres titanium atoms. The 
hexagonal boron net resembles strongly that of graphitic carbon; thus 
we expect strong interatomic bonding within the boron net. The 
titanium atoms nest in interstices provided by the boron net. The 
axial ratio equals 3.228/3.028=1.066, relatively high for the 
AlB 2 -compound group. Each Ti atom is surrounded by twelve equi- 
distant boron atoms. Each boron atom has three boron atoms at a 
short distance, and six titanium atoms at a much longer distance. 



2. Measurements 

2.1 Crystal 

An oriented parallelepiped specimen 1.8 mm x 
2.3 mm x 3.9 mm was prepared from a larger crystal 
grown by a floating zone method [7, 8]. Fluorescent 
x-ray analysis revealed no significant impurities. The 
exact stoichiometry was TiB lt97 . The crystal faces with 
respect to the above dimensions were [2-1-10], 
[-1100], [0001]. These directions represent an a-axis, 
the c-axis, and the direction orthogonal to both. Laue 
x-ray diffraction confirmed these orientations within 1°. 

2.2 Method 

To measure the Cy, we used resonance ultrasound 
spectroscopy, summarized in Fig. 2 [9-11]. Briefly, one 
clamps lightly a regular-shape (cube, cylinder, cube, 
parallelepiped) specimen between two piezo- 
electric transducers. One transducer is swept through 
frequency and the second transducer detects the speci- 
men's macroscopic vibration frequencies (Fig. 3). 
Frequencies of a specimen are determined by five 
factors: (1) shape, (2) size, (3) mass or mass density, 
(4) elastic-stiffness coefficients Cy , and (5) crystal-axis 
orientation relative to macroscopic shape. Thus, by 
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Fig. 2. Schematic measurement setup. Specimen (parallelepiped) is clamped loosely between two piezoelectric 
transducers. One transducer is swept through frequency. The second transducer detects macroscopic resonance 
frequencies, which depend on specimen shape, size, mass, and elastic-stiffness coefficients, the Cy . 
Courtesy of A. Migliori (Los Alamos National Laboratory). 
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Fig. 3. Macroscopic resonance spectrum. Resonance frequnecies/, yield Qj by an inverse-problem calculation. 
Highly overdetermined, the problem uses about one hundred resonance frequencies/, to determine the five inde- 
pendent Cjj. Vertical bars at bottom show predicted resonance frequencies. The inset shows a resonance-peak 
profile, with a Lorentzian shape, the half-power width giving the internal friction Qj\ the imaginary part of the 
total Qj . 



measuring the resonance frequencies f n , one can deter- 
mine by an inverse calculation the stiffnesses Q, . The 
problem is strongly overdetermined: about one hundred 
f n to determine five C y . Not reported here, one can also 
determine the complete internal-friction tensor 
Qi$ l =AfJf n = Qj*/Cij 9 where Af n denotes resonance- 
peak width, Cy* the imaginary part, and Q, the real part 
of the total Q, tensor. Well-described elsewhere [12], 
the inverse problem involves Lagrangean minimiza- 
tion, the Rayleigh-Ritz method, and a least-squares 
procedure for measured and deduced f a values. 

We determined mass density by careful mass and 
size measurements: p = 4.502 ± 0.016 g/cm 3 . This 
compares with reported x-ray mass densities of 4.504 
to 4.53. We ascribe the differences to heavier impurities 
on the Ti sites. From this mass density, we concluded 
our specimen contains no significant void content. 
From the sharp resonance peaks (Fig. 3), we concluded 
our specimen contains few cracks. 



2.3 Errors 

Errors arise from many sources: crystal orientation, 
crystal dimensions, nonparallelism, the inverse prob- 
lem (measured-frequency to C ijkl conversion), speci- 
men-transducer interactions. Several authors described 
these errors elsewhere, especially the first reference 
[13-15]. The effect of the slight (1 %) departure from 
stoichiometry is hard to estimate. Because of the strong 
covalent bonding within the boron planes, we conjec- 
ture a very small error arising from an occasional miss- 
ing boron atom. 



3. Results and Discussion 

Table 1 shows our principal results: the Q, and their 
uncertainties. (Because the two previous reports omit- 
ted error estimates, very detailed comparisons are 
precluded.) Only five Q, are independent because 
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Table 1. Monocrystal and polycrystal elastic constants of titanium 
diboride. Unless specified, all units are GPa, except for the Poisson 
ratios, v, which are dimensionless 



Table 1 also shows the principal Young moduli E i{ 
computed by 





Present 


Gilman - Roberts 


Spoor et al. 


p(g/cm 3 ) 4.502 ±0.016 








Monocrystal elastic constants 








Cn 


654.4 ±1.9 




690 


660 


C33 


458.1 ±1.4 




440 


432 


C44 


262.6 ± 0.3 




250 


260 


Qe 


302.7 ±0.7 




140 


306 


Q2 


48.98 ±1.4 




410 


48 


C13 


95.25 ±0.55 




320 


93 


En 


633.3 ±2.2 




389 


639 


^33 


432.3 ±1.7 




254 


408 


v 12 


0.0460 ±0.0022 




0.3877 


0.0437 


v 13 


0.1984 ±0.0015 




0.4453 


0.2059 


^31 


0.1354 ±0.0010 




0.2909 


0.1314 


Debye 


characteristic temperature 







D (K) 1217 ±6 



1211 



Voigt -Reuss-Hill- average quasiiso tropic (polycrystal) elastic constants 
C L 599.9 642 593 

B 247.5 417 244 

G 264.3 169 262 

E 584.7 446 579 

v 0.1063 0.3219 0.1037 



a Ref.5. 
b Ref. 6. 



Q6 = (Qi ~ C^/2. As usual (they fail to correspond 
directly to a phonon), the off-diagonal C y - (C 12 and C 13 ) 
show the largest uncertainties. Table 1 also gives the 
Spoor et al. results [6], which differ from ours by an 
average of 2.2 %. The average uncertainty in our six C y - 
is 0.7 %. The close agreement between our results and 
the Spoor et al. results suggests strongly that the earli- 
er Gilman-Roberts results [5] are wrong, the largest 
discrepancies occurring in the off-diagonal Q, : C 12 and 
C 13 . Further support for the correctness of the 
Ledbetter-Tanaka/ Spoor et al. results arises from 
ab initio calculations that yielded 5 = 251 GPa and 
v= 0.12 [16], versus 417 and 0.32 for Gilman-Roberts. 



E *=% 



0) 



Here S^ denotes the elastic-compliance tensor, the ten- 
sor inverse of the C y - . As we expect from C n > C 33 , E n 
exceeds E 33 ; that is, TiB 2 is much stiffer within 
the basal plane than along the c-axis. Obviously, this 
relates to the crystal structure where the covalent- 
bonded boron atoms lie in the plane perpendicular 
tox 3 . 

Table 1 shows also the three principal Poisson ratios 
Vy computed by 



-s./s. 



(2) 



Within the boron plane, the Poisson ratio v 12 is 
extremely low, reflecting strong covalent bonding and 
the strong resistance of boron atoms to change their 
bond angles. The v 13 Poisson ratio is only slightly 
below normal, indicating weaker bonds out of the 
boron-atom planes than those within the planes. 

The shear elastic anisotropy of hexagonal crystals 
can be expressed in various ways. The simplest is 
C 66 /C44, 1.15 for TiB 2 , thus weak elastic anisotropy 
(the isotropic case corresponding to 1.00). Because the 
Young modulus depends so strongly on the shear mod- 
ulus, a Young-modulus variation with direction also 
serves as an effective shear-anisotropy indicator. Spoor 
and colleagues showed a Young-modulus representa- 
tion surface; it is nearly spherical [6]. From the alter- 
nating-layer boron-titanium crystal structure (Fig. 1), 
one expects higher elastic anisotropy than one finds. 

The lower part of Table 1 gives the averaged- 
over-direction quasiisotropic elastic constants obtained 
from the C y by a Voigt-Reuss-Hill average [17]. These 
constants include longitudinal modulus C L , bulk mod- 
ulus B, shear modulus G, Young modulus E, and 
Poisson ratio v. These are the elastic constants appro- 
priate to a full-density nontextured polycrystalline 
aggregate in which the grain boundaries cause no soft- 
ening. Grain size produces no effect on elastic con- 
stants if it is small relative to specimen size. (A review- 
er pointed out that nanosize grains may soften the 
elastic constants.) The most unusual feature of our 
results is the high GIB ratio, thus low Poisson ratio. 
Our nonpublished results on a YB 66 monocrystal gave 
an averaged-over-direction (Kroner method) Poisson 
ratio of 0.13. For boron, a handbook value is v= 0.089 
[18]. Both the bulk and shear moduli of TiB 2 
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exceed considerably (by 30 % to 40 %) the handbook 
values of boron: 248/179 = 1.39 and 264/203 = 1.30, 
showing the strong interatomic bonds in TiB 2 , both for 
extension-compression and for shear-bending. Beside 
implying a low Poisson ratio, the very high GIB ratio 
holds implications for several other crystal-bonding 
properties: covalency, Cauchy-relation departure, 
many-body forces, and others. We plan to discuss all 
these elsewhere. 

The elastic stiffnesses yield three useful sound 
velocities: the longitudinal velocity 



v 1= (C L /p) 1/2 = 



S + (4/3)G 



(3) 



the shear or transverse velocity 



and the mean velocity (as defined by Debye) 



_3___2_ J_ 



(4) 



(5) 



We found v^l.54, v s = 0.766, and v m = 0.835 cm/jLis. 
From the C§ and the mass density, we can compute 
the acoustic Debye characterisitic temperature <9 D . At 
zero temperature, the acoustic <9 D becomes identical 
with the calorimetric 6^ [19]. O u is proportional to the 
mean sound velocity: 



0n=Kv m , 



(6) 



where v m denotes mean sound velocity and K is given 
by 



,( 



K = - 

k 



r 



AnV 



(7) 



Here h denotes Planck's constant, k Boltzmann's 
constant, and F a atomic volume. The velocity v m comes 
from the integration over all directions: 



^t=jl^dQ/4K. 



(8) 



a=l,3 



Here v 1 denotes the quasilongitudinal wave velocity, 
v 2 and v 3 the quasitransverse wave velocities, and d£l 



an increment of solid angle. Equation (8) can not be 
integrated analytically, and numerous numerical and 
approximate methods have been used for its solution. 
Phase velocities v a are roots of the Christoffel 
equations: 



det(C ijkl x j ^-pv 2 5 il )=0. 
(sum over repeated indices) 



(9) 



This expression follows from equations of motion 
for plane, monochromatic waves, where p denotes 
mass density, C ijkl fourth-rank elastic-stiffness tensor, 
jCi components of the unit wave vector relative to the 
cubic axes, and 8- A the Kronecker operator. Equation (9) 
usually yields three distinct real roots pv 2 . 

We computed an exact <9 D by using a distribution of 
489 vectors proposed by Overton and Schuch [20], 
which we distributed over the usual 48 [100]-[1 10]- 
[111] stereographic unit triangles, thus a total of 23 472 
directions. We obtained <9 D =1217K. This result 
differs enormously from the handbook calorimetric 
value, 1576 K. Often, calorimetric values contain large 
errors because of large extrapolations to zero tempera- 
ture and/or large uncertainties in the lowest-tempera- 
ture lattice specific heat. 

From B and from handbook values of heat capacity 
C, volume thermal expansivity j8, and volume F, we 
computed the effective Gruneisen parameter: 



y = B s pV/C r 



(10) 



We obtained y= 1.71, the handbook values for boron 
and titanium being 1.85 and 1.33. Alternative gammas 
can be computed when one knows the third-order 
elastic-stiffness coefficients, the C ijklmn [21]. The few 
gammas known for compounds similar to TiB 2 pre- 
clude any comparisons. Following Pearson's reasoning 
[22], the Ti-B bonds would lead to y=2. 

Finally, because some authors suspect voids/cracks 
in TiB 2 , we want to describe briefly how the above full- 
dense quasiisotropic elastic constants would change 
with voids. Focusing on A1 2 3 , Ledbetter, Lei, and 
Datta [23] gave a theory for void effects on elastic 
constants. Principal results include the following: 
Voids soften the bulk modulus more than the shear 
modulus. In the dilute limit, for spherical voids, 
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our results agree with the classical results of Mackenzie 

[24]: 



M = _ 3C; (v =l/3) 



&G 15 

- = - 7 c. (v =l/3) 



(11) 



(12) 



Here, c denotes void volume fraction and v the 
void-free-state Poisson ratio. However, we emphasize 
that the often-used rule-of-thumb that elastic stiffness 
varies as mass density is true only in the dilute limit. 
Void shape plays a key role, especially if the voids 
possess an oblate-spheroid (disc) shape. Dunn and 
Ledbetter [25] focused on the interesting, unexpected 
effects of voids and cracks on the Poisson ratio. Other 
authors addressed this problem using various approxi- 
mations [26]. 



4. Conclusions 

1. Our elastic-stiffhess-coefficient measurements on 
titanium diboride support the report of Spoor and 
colleagues rather than the older (now handbook) 
values of Gilman and Roberts. The principal differ- 
ences in the two previous Q,- sets lie in C n and C 13 . 

2. From the Cy we computed several additional use- 
ful physical properties: 

(i) sound velocities; (ii) Debye temperature; 
(iii) Gruneisen parameter. 

3. Our computed acoustic Debye temperature, 1217 K, 
is about 20 % lower than the handbook calorimetric 
value. 

4. Our computed Gruneisen parameter, 1.71, suggests 
the importance of Ti-B bonds along with B-B 
bonds. 

5. Using the C^/C^ ratio as a shear-mode elastic- 
anisotropy criterion, TiB 2 shows low shear- 
mode anisotropy 1.15. A surprise because of the 
B-Ti-B . . . layer crystal structure. 

6. Our elastic constants estimated for a full-dense 
polycrystal depart strongly from those given in 
Munro's review [2], especially in the uncertainties. 
For example, Munro proposed a 70 % uncertainty 



in the Poisson ratio and a 24 % uncertainty in the 
bulk modulus. Our polycrystal results differ from 
the Spoor et al. result by only 2 %. This finding 
agrees with the well-known fact that modern 
measurement methods give the elastic-stiffness 
coefficients easily within one percent. 

7. Finally, we are surprised by the low elastic shear 
anisotropy shown by such a strongly layered 
crystal structure. Reference [8] gives a possible 
explanation. Part of the boron-atom />-obitals lie in 
a p-d hybridized band, weakening the j^-electron 
contribution to B-B bonding. If the bonding weak- 
ens, the anisotropy decreases. 
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